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Abstract 

We prove the existence of Absolutely Minimizing Lipschitz Extensions by 
a method which differs from those used by G. Aronsson in general metrically 
convex compact metric spaces and R. Jensen in Euclidean spaces. Assuming 
Jensen's hypotheses, our method yields numerical schemes for computing, in 
euclidean R'\ the solution of viscosity of equation Aoo(m) ~ with Dirichlet's 
condition. 



1 Introduction 

To produce an optimal solution to Tietze's extension problem in general metrically 
convex compact metric space, we have introduced a class Hh of extension schemes 
which solve the problem |7]. 

In this paper we first prove in section 3 that, for any continuous Dirichlet's condi- 
tion /, there exists a subsequence (^^/i(n)(/))nGN which converges to an AMLE of 
/. Therefore, assuming Jensen's hypotheses 0], H^if) approachs the solution of 
viscosity of Aoo{u) = under Dirichlet's condition / when h tends to 0. 

Unfortunately it is generally hopeless to obtain a numerical approximation of this 
solution on, say, a regular grid of step h by discretisation of Hh on this grid. In 
fact, by such a discretisation we obtain an extension which is Lipschitz-optimal 
not for euclidean metric but only for the geodesic metric on the grid. 

To overcome this difficulty we introduce in this paper an explicit scheme of ex- 
tension valid on any finite network contained in the considered metrically con- 
vex compact metric space and prove that the extension converges to an AMLE 
when the network suitably densifies the metric space. As a consequence, assuming 
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Jensen's hypotheses j^, we obtain numerical approximations of the solution of vis- 
cosity of Aoo('u) = under Dirichlet's condition. Note here that A. Oberman [S| 
has obtained very similar numerical approximations in R", based upon the same 
numerical scheme, proposing a proof of the convergence of the scheme based upon 
the Aoo— approach of the problem. 

In the whole paper {E, d) denotes a metrically convex compact metric space that 
is a compact lenght space with the terminology of ([2], appendix). We denote by 
S the Hausdorff metric induced by d on compact non-empty subsets of E. 
The second part of the paper is organized as follows. 

In section 4 we prove that solutions of Hl.lfl (see below) satisfy the maximum 
principle and, as a corollary, uniqueness of the solution. 

In section 5 we prove the existence of the solution of 1)1. 1|1 and we study the 
stability of this solution. 

In section 6 we prove the existence of an AMLE as the limit of solutions of l|l.l|l 
for sequences ((G^, Ki))neN which suitably densify E. 

Definition 1.1. A network on i? is a couple (G, V) where G C E denotes a finite 
non-empty subset of E and V a mapping x G G ^ ^i^) C G, {V{x) is the neigh- 
bourhood of x) which satisfies 
(PI) for any x € G, x € V{x); 
(P2) for any x,y e G, x € V{y) iff y G V{x); 

(P3) for any x,y £ G, there exists xi,X2, Xn-i, Xn € G such that xi — x, Xn — y 
and Xi € V{xijf-i) for i = 1, n — 1; 

(P4) for any a; G G, any y e G — V{x) there exists z G V{x) such that d{z,y) < 
d{x,y). 

To any chain such as in (P3) we associate its lenght d{xi, Xi+i). We define 

the geodesic metric dg on (G, V) by letting dg(x,y) be the infimum of the lenght 
of chains connecting x and y. 

It follows from (P1),(P2),(P3) that dg is a metric, that d{x,y) < dg{x,y) for 
x,y G G and that d{x, y) — dg{x, y) for x,y £ G, x € V{y). It follows from (P2),(P3) 
that if G has at least two elements (assumed from now on) then V{x) — {x} ^ 
for any x £ G. We shall denote V{x) :— V{x) — {x}. Extra-condition (P4), crucial 
in this paper (see the end of theorem 14.11 and theorem 16.31 iii) , will be used as 
follows: 

for any x £ G, D non-empty subset of G, d{x, D) > 0, there exists y £ V{x) such 
that d{y,D) < d{x,D). 

We consider the following functional equation with Dirichlet's condition : 



J u{x) = iJ,{u; x) yx £ G — S: 
\ u{s) = f{s) Vs G S. 



(1.1) 
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Here S denotes a non-empty subset of G, function / is the Dirichlet's condition 
defined on S, u is the numerical unknow function defined on G and 

fj.{u;x) = inf sup M{u;z,q){x) ; (1-2) 

zeV(x) q^V{x) 

where 

d{x,z)u{q) +d{x,q)u{z) 

■= d(.,z) + rf(x,g) • ^'-^^ 

Remark 1.2. It can be checked that 

IJl(u;x) = sup inf M(m;z,(7)(x) . (1-4) 

2ev(2;)9e^(^) 



It can also be checked that 



J{li{u]x)) = inf J(m) 



where 

T/ / I - ^^ I 

J(/x(u;a;))= sup — — . 

zev(x) (J-\x,z} 

Therefore ^(u; x) is the explicit solution of the problem of minimization considered 
by A.Oberman. 



2 Basics 

Let / be any function from dorn{f) C -B to M. We define by 
K{f) := sup — — . 

X ,y^dom{f) ,x^y CLl^X^ y) 

We call concave modulus of continuity any mapping oj : M+ IR.+ which satisfies 
the following: 

(i) cij(O) = and lu is continuous at 0; 

(ii) uj is increasing: hi < h2 < w(/i2); 

(iii) u! is concave. 

We say that / is i7— continuous iff there exists a concave modulus of continuity uj 
such that, for any x,y ^ dom{f), 

\fix)-f{y)\<u;{dix,y)). (2.1) 

For such a function /, we denote by i-u{f) the lower bound of those concave moduli 
of continuity which satisfy (|2.1|) . 

For any A C dom{f), we have obviously u!{f \ A) < uj{f) (symbol | denotes 
restriction to). 

Let us restate here results of [S] which are of constant use in this paper. 
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Proposition 2.1. Let f, g be any two ^-continuous real-valued functions of do- 
main S and let A and B he any two compact non-empty subsets of S. Then 

ll^(/)-w(.9)lloo.K+<2||/-<?||oo.s; (2.2) 

\\ u:{f \ A) - u{f \ B) |U,R+<4c.(/;<5(AS)). (2.3) 

Note that H2.2|l and p. 3(1 have been estabhshed in |Bj for weak moduH of continuity. 
It is immediate that these inequahties hold for concave moduh of continuity with 
the same constants. 

Remark 2.2. So, aside the obvious fact that AMVlE (see below) are more general 
than AMLE, the true reason why we adopt the modulus of continuity approach 
rather than the Lipschitz approach in this paper is that there is no equivalent of 
(|2.2|l and (|2.3|l for Lipschitz functions. 

Now we recall Aronsson's definition of an AMLE jj. Let e be a Lipschitz extension 
of a Lipschitz function / of compact domain. 

Definition 2.3. We say that e is an Absolutely Minimizing Lipschitz Extension 
of / if for every non empty open D C E, D D dom{f) — $ we have 

K{e I D) = K{e I dD), 

where dD denotes the boundary of D. 

Characterisation below has been noticed by Aronsson Pp. 

Proposition 2.4. An extension e of f is AMLE iff for any non empty open 
D d E, D n dom{f) = 9 we have 

e{x) < inf {e{y) + n{e \ dD)d{x,y)) ,yx e D, (2.4) 

yedD 

and 

sup {e{y) - K{e \ dD)d{x, y)) < e{x),Wx G D. (2.5) 
vedD 

In this paper we use a slightly more general definition (see remark . Let e be 
a continuous extension of a fi-continuous function /. 

Definition 2.5. We say that e is an Absolutely Minimizing fl Extension of / if 
for every non empty open subset D of E, D n dom(f) — we have 

w(e I D) = uj{e \ dD), 



where dD denotes the boundary of D. 
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The analog of proposition 12 . 41 is : 

Proposition 2.6. An extension e of f is AMflE iff for any non empty open 
D d E, D n dom{f) = we have 

e{x) < inf {e{y) + Lu{e \ dD; d{x, y))),yx £ D, (2.6) 
t/eS-D 

and 

sup {e{y) - Lu{e \ dD; d{x, y))) < e{x),yx e D. (2.7) 

yedD 

It follows from these definitions that, if / is Lipschitz and e is an AM^IE of /, 
then e is an AMLE of /. 



3 Convergence of harmonious extensions to an 

AMQE 

Let / be a continuous function of closed domain dom(f) C E. 

For any h > 0, for any x £ E we denote by Vh{x) the closed ball of center x, radius 

r{x) = inf(/i, d{x, dom(f))). 

By the [theorem 3.3 of [Jj], there exists a unique continuous extension Hh{f) from 
£^ to R which satisfies functional equation 

5(x) = i sup g{z) + ]- ini g{z),yx£E. (3.1) 

^ zeVh(x) z zeVh(.x) 

Moreover, as noticed in remark 3.4 of [J], the proof of [theorem 3.3 of ^] shows 
that 

c.(ff,(/))-c.(/). 

Lemma [3. II below shows that Hh{f) is close to an AAIilE of /. Its proof uses the 
arguments of Proposition 3.9 of 

Lemma 3.1. For any non empty open subset D of E , DC]dom{f) = 0, we have 
Hh{f){x) < inf {Hh{f){y)+uj{Hh{f) \ dD;d{x,y))) + 2w(/;/i),Vx G D, (3.2) 

yEau 

and 

sup {Hh{f){y) - u;{H,{f) \ dD; d{x, y))) - 2u;{f; h) < Hh{f)(x),yx G D. (3.3) 

yedD 

Proof. Since the arguments are symmetric we prove only H3.2|l . 

Since i5 is a compact metrically convex metric space, by theorem 3.3 'T there 
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exists a unique extension v of Hh (/) | dD in E such that 

v{x) — — sup v{z) + — inf v{z), Vx G i?. (3-4) 

where ^^^(x) := {z e £^ : d{x,z) < inf{h,d{x,dD))}. 
Moreover 

- viy) < LoiHhif) \ dD; y)), \/x, y e E. 
In particular we have 

v{x) - v(y) < uj{Hh{f) I dD; d{x, y)), Vy £ dD, Vx E D. (3.5) 

Now, let us bound sup^g^, | Hh{f){x) — v{x) \. By symmetry we have only to 
bound from above: A — sup^^jj{Hh{f){x) — v{x)). Let 

F^{xeD: HhiDix) - v{x) = A}, M - sup Hh{f ){x), 

and 

F = {x e F : Hh{f){x) = M}. 

Let a; G be such that 

d{x,dD) = inf d{y,dD). (3.6) 

yeP 

Let us first show that we cannot have d{x, dD) > h. Towards a contradiction let 
us assume it is the case. 

Then we have Whix) — {z £ E : d{x, z) < h}. 

Since dom{f) Ci D — ^ we infer that d{x, S) > d{x, dD) > h, that is Vh{x) = {z e 
E : d{x,z) < h}. Therefore Vh{x) = Whix). 

Now, using a similar argument to this of theorem 3.3-uniqueness- of [7], we infer 
that Vh{x) C F. It follows that there exists z € F such that d(z,dD) < d{x,dD) 
which is a contradiction with definition (|3.t)|l of x. 

Now, since d{x, dD) < h, there exists y G dD such that d{x, y) < h and Hh{f){y) — 
v{y)- 

By f2-stability of both v and Hh{f) we have: 

Hhif){x) - v{x) = Hh{f){x) - Hh{f){y) - v{y) - v{x) < 2io{f; /i).Therefore 

A<2t^(/;M- (3.7) 

Now \ei X E D and y G dD. We have 

Hh{f)[x) - Hh{f){y) ~ LuiHhif) I dD; d{x, y)) = A, + 

where 

A^^Hu{f){x)-v(x), 
and 

A2 = v{x) - v{y) - uo{Hh{f) \ dD; d{x, y)). 

From inequahties ifT^I) and iflTTjl . we obtain Al + A2 < 2uo{f; h) + 0. □ 
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Now we prove the existence of AMVIE in any metrically convex compact metric 
space. 

Theorem 3.2. Let {h{n))ni=n be a sequence of positive reals which converges to 
0. If sequence (-ff/i(n) (/))neN converges uniformly to a continuous extension g of 
f then g is an AMflE of f . 



Proof. By symmetry we prove only H2.6(l . 

Let D be non empty open subset of E such that D D dom{f) ^ 0. For any e > 0, 
there exists > such that Vn > N , \\ H^n){f) ^ 9 ||oo.£;< £• For any x ^ D, 
y £ dD, n > N we have 

g{x) - g{y) - uj{g \ dD; d{x, y)) < Ai + A2 + A3 + A^, where 
Ai =1 g{x) ~ Hh(^^){f){x) |< e; 

A2 =1 Hh(n){f)(x) - Hh(n){f){y) - ^{Hh(n){f) I dD;d{x,y)) |< 2uj{f;h{n)); 
A3 =1 u;{Huin){f) I dD;d{x,y))~uj{g \ dD;d{x,y)) \< 2e; 

^\ H^in)if)iy) - giy) l<e. 

The second inequality follows from Lemma f3. II and the third one from (|2.2|l . We 
obtain g{x) - g{y) - uj{g \ dD] d{x, y)) < 4e + 2w(/; h{n)). 
By letting n — > 00 we have 

gix) - g{y) - uj{g \ dD; d{x, y)) < 46, Wy G dD, 

and we obtain the stated result. □ 

Theorem 3.3. For any continuous real-valued function f whose domain is a 
compact non-empty subset of E, there exists an AMflE of f . 



Proof. The set {Hh{f),h > 0} is equicontinuous and equibounded. Therefore, 
by Ascoli's theorem, there exists a subsequence (i?/i(n)(/))neN which converges 
uniformly to a continuous extension of / which is an AMnE of / by theorem 
El □ 

Remark 3.4. Moreover if, for any /, there exists a unique AMilE of / denoted 
by H{f) then limh^o Hh{f) = H{f). In this case it follows from proposition 3.9 
of 13 that 

II H{f I A) - H{f I B) U.E< Mf: S{A, B)), 
for any non-empty compact subsets A,B of dom{f). 

Remark 3.5. We can summarize the difference between Jensen's proof (4j of the 
existence of an AMLE and our own proof as follows. Jensen obtains the desired 
AMLE as a limit of local (because solutions of PDE) extensions which become 
more and more optimally Lipschitz. We obtain the desired AMLE as a limit of 
optimally Lipschitz extensions which become more and more local. 
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Aronsson (see also |S] and [S] ) proves the existence of AMLE by giving two 
explicit solutions: 

u = sup{w : w AMLE of / from above }, 

V = inf{w : w AMLE of / from below }. 

Our proof leads to less explicit but, assuming uniqueness, more constructive solu- 
tions than Aronsson's one. 

Remark 3.6. Note the formal analogy between the process u — > ^h{u) of harmo- 
nious regularization defined by 

<i>h{u){x) sup u{y)+ inf u{y)) 

which deals with PDF AoqU — and the process u of harmonic regu- 

larization defined by 

^h{u){x) = — ^ — 

which deals with PDE Au^O. 

It is known since Gauss that any harmonic function satisfies — u for any h > 

0. The analog of this result does not hold in general for the process of harmonious 
regularization : it can be seen by numerical tests that <i>/i(u) ^ u for u{x,y) ~ 
x'^^^ — y'^^^ even in subdomains where this function is analytic. 

However some functions u solutions of A^oU = have this property : for example 
linear functions, {x\ -I- a;|)^/^, arctan{xi/x2), in euclidean plane, xf — x\ in the 
plane equipped with sup norm. 

Remark 3.7. When {E,d) is || II2) with open convex non empty subset 
of euclidean M", dom{f) = dfl, it can be shown directly (that is whithout us- 
ing theorem 13.31 and the equivalence between AMLE and solution of viscosity of 
AooU = 0) that Hh{f) converges, when h tends to 0, to the solution of viscosity 
of AooU = 0, u I dfl = f. It is a consequence of Jensen's uniqueness results ^ and 
of a Barles-Souganidis's result [3] : see Appendix. 

Remark 3.8. The results of fTI and of this section hold for spaces more general 
than compact metrically convex metric spaces. They hold in compact metric spaces 
{E, d) having the following properties: 

i) 

- sup inf d{q,r) + — sup inf d{r,q) < d{x,y), for any x,y £ E 

'2 qeBhix)r&Bh.{y) r£B,,(y) 1^Bh(x) 



ii) For any x € E and y ^ Bh(x) there exist z £ Bh{x) such that d{y, z) < d{y, x). 
It follows that the condition of convexity on fl assumed in remark 13.71 can be 
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removed. 



Note that conditions i) and ii) can hold in metric spaces which can be very far 
from metrically convex metric spaces (some finite metric spaces satisfy conditions 
i) and ii)) : we have therefore established a theorem of existence of an AMLE 
under weaker hypotheses than those of Milman ^ and Juutinen (however our 
result holds only for compact spaces). 

4 Uniqueness theorem for functional equation ( 11.11) 

As usual, we first prove a maximum principle. 

Theorem 4.1. Let f , g he any two real-valued functions both of domain S . Let 
u,v he two solutions of 1^1.1]) with Dirichlet's conditions f and g respectively. Then 



Proof. Let us set A = sup^^q{u{x) — v{x)), F ~ {x £ G : u{x) — v{x) = A}, 
A = supj-g^ u{x) and F = {x d F : u(x) — A}. 



sup(u(a;) - v{x)) < sup(/(s) - g{s)). 



(4.1) 




We start our proof by choosing x G F such that d(x, S) — inf ^ 
If d{x, S) = then equality H4.1|l is true. 



yeF 



d{y,S). 



Else, let zi e V{x) such that 



sup M{v; zi,q){x) = n{v;x) . 



We have 



A< sup M{u; zi,q){x) — sup M{v;zi,q){x) . 

qeV{x) qeV{x) 



Let qi e V{x) such that 



M(u;zi,gi) 



sup M{u; ^1, (?) . 



We have 



A < 



d{x, zi){u{qi) - v{qi)) ^ d{x, qi){u{zi) - v{zi)) 
d{x, zi) + d{x,qi) d{x, zi) + d{x,qi) 



< A 



from which it follows that 2:1,(71 G F. 
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Since 

u{x) = ^{u;x) < sup M{u; zi,q){x) — M{u; zi,qi) < u{qi), 
qev{x) 

we have qi € F and u{qi) = u{x). 
Since 

- M TTTTT ^ ' 

a(x, zi) + d[x, qi) 

we have u{zi) = u{x). 
Since 

u{x) — u{qi) — sup M{u;zi,q) 

and 

w ^ t7/ n / \ d{x,zi)u{q) +d{x,q)u{zi) 

Vq e , u(a;) > M{u; zi,q) = — — ^ 

d[x, zi) + d(x, q) 

with u{zi) = u{x), 
we have 

Vz G V{x) , u{z) < u{x). 

We finish the proof of our assertion by prooving that u is constant in V{ 
Towards a contradiction let q G V{x) such that u{q) < u{x). We have 

d{x,t)u{q) + d{x,q)u(t) 

u{x)< sup M(?/;g,t)= sup ^ . • 

tevix) tGV{x) a[x,i) + a(x,q) 

Using (|4.3|) . we obtain 

. ^ ^ dfa;, + g)u(a;) 
< sup . 

tev{x) d{x,t) + d{x,q) 

Let c > such that u{q) — u{x) — c, we have 

dfa;, — c) + (i(a;, (7)u(a;) 
u[x) < sup ,^ , ^ • 

tevix) d{x,t) +d{x,q) 

On the other hand, we can write 

x) < sup {u{x) - -) , 

t(^v{x) d{x,t) + d{x,q) 

that is 

inf <0. 

tev{x) d{x,t) +d{x,q) 

Since 

d{x,t) 

mf — J- > , and c> 0, 

t(^v{x) d{x,t) +d{x,q) 
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we obtain the desired contradiction. 
Now, using inequality 

A< sup M{u;zi,q)~ sup M{v;zi,q), 

qeV(x) q£V{x) 

and u = constant in V{x), we have 

d{x,zi){u{q) - v{q)) d{x, q){u{zi) - vizi)) 



\/q G V{x),A < 



d{x , zi) + d{x , q) d{x , zi) + d{x , q) 



Therefore A < u{q) - v{q) < A, Vg e V{x), that is V{x) C F. 

Since u is constant on V{x) we have V{x) C F. Since V satisfies (P4), we have 

d{x, S) = iniy^p d{y, S) and M^^y^^^ d{q, S) < d{x, S) 

which is a contradiction with V{x) C F. So (|4.1|) is proved. □ 

As immediate consequences of theorem 14.11 we obtain Theorem 14 . 21 and coroUary 
Ol 

Theorem 4.2. Functional equation 11 .1]) has a unique solution. 
Corollary 4.3. Let u a solution of fi. j|) then 

inf /(s) < u{x) < sup f{s) ,\/x £ G. (4.4) 

5 Existence and stability of solutions of ( 11.11 ) . 

To prove the existence of a solution of (|l.l|l . we introduce a process of evolution 
u — > whose the stationary state u — is solution of (|l.ll) . Precisely 

= \E'(m; a; at) o \I>(m; sat^i) o ... o \I>(u; a;i) where {xi, xat} is an enumeration 
of G — 5 and ^'(u; x), x £ G — S is defined as follows: 

■^{u;x){y) =u{y) \iyeG-{x}; 
^[u]x){x) — ^{u;x) ii y = x 

We need three lemmas useful for existence and stability. 

Lemma 5.1. For any two scalar-valued functions u,v of domain G, we have the 
following properties : 

u<v =^^{u]x) <^{v;x) ,'ix £G - S ; (5.2) 



I ^{u- x){y) - -^{u- x){z) |< u{u- dg{y, z)) , \/x e G - S ,\/y, z^G ; (5.3) 
sup(*(m; x){y) - *(u; x){y)) < sup{u{y) - v{y)) . (5.4) 

y6G ySG 
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Proof. Let us show (|5.2|) . Let u, v scalar- valued functions of domain G such that 
u < V. Let X € G — S. We have 

^{u\x){y) = u{y) < v{y) = for e G - {x}. 

Since Vz, 5 G V^(a;) we have M(u\z.,q)(x) < M{v] z,q){x), therefore fi{u;x) < 
IJ,(v;x) 

Let us show (|5.3|l . 

It suffices to prove that | ^{u;x){x) — u{y) |< uj{u; dg{x,y)) , Vy G G. 
Let y e G — {x} we have two case : 

First case : suppose that y e G — V{x). 

Let zi € V{x) such that dg{x, y) — dg(x, zi) + dg{zi, y). 

We have 

^[u,x){x)-u{y)< sup + -..(y)) ^ 

By definition of uj{u) we have 

*(.;x)(x)-^y)< ( d{x, z,Mu; dgjg, y)) + d{x qMu; dgjz^, y))) 

q(zv(x) d{x,zi) +d{x,q) 

By concavity of uj{u) we have 

,T,/ ^/^ ^ d{x, zi)dg{q,y) + d{x,q)dg{zi,y) 

^'(u;x)(a;) - < sup f — — ), 

gev{x) d{x,zi) + d{x,q) 

Since 

d{x, zi)dgd{q, y) + d{x, q)dg (21 , y) = zi ) {dg (q, y) ~d{x,q))+ d{x, q)dg (x, y) 
by the triangle inequality, we have 



d{x,zi){dg{d{q,y)) + d{x,q){dg{zi,y))) 
d{x, zi) + d{x, q) 

Second case. Suppose that y G V{x). We have 



< dg{x,y) . 



^{u;x){x)-u{y)< sup {M{u;y, q) - u{y)) < sup dixyMu,diq,y)) 
qev(x) qev(x) d{x,y) + d{x,q) 

By concavity we have 

nu;x){x)-u{y)< sup ^(^. y^.^Kfa;/) ). 

q€Vix) dg{x,y)+dg{x,q) 
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Since 

dgil, y) ^ ^ 

dg{x,y) +dg{x,q) ~ 

we conclude that "i'{u;x){x) — u{y) < uj{u;dg{d{x,y))) . The arguments to prove 
that u{y) — ^{u]x){x) < uj{u; dg{d{x,y))) are symmetric using (|1.4|l instead of 
l|1.2|l . Inequahty (|5.3|l is therefore proved. 

Inequahty H5.4|l holds because we have x) — fi(v; x) < sup^gy(-^-j(M(z) — v{z)). 

□ 

Lemma 5.2. For any scalar-valued functions u,v of domain G, we have the fol- 
lowing properties : 

u<v=^ < ,\fxeG- S ; (5.5) 
I ^u){y) - $(w)(z) |< uiu; dg{y, z)) , yx e G ^ S , \/y, z e G ; (5.6) 



and 



inf u{z) < ^{u){x) < supu(z) , Va; G G . (5.7) 



Proof. The proof is a consequence of Lemma l5. II □ 

Now let Uo be defined by 

Uo{x) = m{{f{s)+Lu{f;dg{x,s))) , Vx £ G . 

Function Uq looks like classical Af "^Shane maximal Lipschitz-optimal extension of 
/ on G. But here Uq is defined with both d (in w(/)) and dg {in dg{x, s)). Therefore 
we have to check that Uq is an extension of /. 

Lemma 5.3. We have 

Uois)^fis),forseS; (5.8) 

I Uoix) - Uo{y) \< u{f-dg{x,y)), for x,y e G ; (5.9) 
inf /(s) < Uo{x) < sup/(s), forxeG . (5.10) 



Proof. Let s G we have 

UoiS) - fi~s) < fiS) - Loif; dgi~s, ~s)) - f{S) < ; 
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and 

/(s) - Uoirs) = sup(/(S) - f{s) - io{f; s))) , 

therefore 

- UoiS) < sup(L^(/; s)) - ^(/; s))) . 

Since d < dg we have /(s) — ?7o(s) < and /(S) = Uo{s). 
Let a;,y G G we have 

Uoix) - Uoiy) < sup(/(s) + ^^(a;, s)) - /(s) - dg{y, s))) , 

s£S 

Therefore 

Uo{x) - t/o(y) < sup(w(/; dg{x, s)) - dg(y, s))) < dg(a;, y)) . 
ses 

Let X e G we have 

L/o(x) < inf /(s) + suptj(/;dg(a;,s)) . 

since a)(/;dg(a;,s)) < sup^^_^^gs.(/(si) - /(sa)), we have 

C/o(a:) < inf /(s) + sup (/(si) ~ /(sa)) < sup/(s). 

^£•5 si,s2es s£S 

Since (ig(a;, s) > 0, we have Uo{x) > inf^gg /(s)- 

□ 

Now we are ready to prove the existence of a sohition of l|l.l(l . 

Theorem 5.4. Let {Un)neN the sequence defined inductively by Un+i = $(CX,i), Vn G 
N. T/iis sequence converges to a solution of denoted by K{f): 

K{f){s)^f{s),yseS (5.11) 



Moreover we have 

I K{f){x) - KU){y) |< c.(/; dg{x, y)) , ^x,y e G 



(5.12) 
(5.13) 
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Proof. Let us show that {Un)nGN is decreasing. By Lemma it is sufficient to 
prove that Ui < Uq. Given an arbitrary k G {1, , N}, we have 

'^{Uo;xk)ixk) - Uaixk) ii{Uo\Xk) ~ Uo{xk). 

Let Sk & S such that 

/(sfc) + uj{f; dg{xk, Sk)) = inf (/(s) + uj(f; dg{xk, s))). 



Since 



inf (/(s) + < f{sk)+uj{g){dg{y,Sk)) , Vy £ V{xk) 



we have 

nuo;xk)ixk)-fisk)< inf sup ( ^.(^^-^ ^)^(/; ^fe)) + y-^.-. ^.(^. ^.)) . 

zGV(i;fc)gg(>(^^) dg{Xk,z}+dg{Xk,q) 

By concavity we have 

ir.frr \( \ fi \ ^ ■( ,rdg{xk,z)dg{q,Sk) + dg{xk,q)dg{z,Sk). 
'i>{Uo;Xk){xk)~f{sk) < mf sup cj(/; , /, ^-^^ ). 

First case. If Sfe ^ V"(a;fc), then 3z' e V^(a;fe) such that dg{xk,Sk) — dg{xk,z') + 

dg{z',Sk). 

We have 

nuo;xk)ixk)-fisk)< sup ^ 

ggy(^^) dg(xfc,2;') + dg(Xfe,g) 

Since 

dg{xk,z')dg{q, Sk)+dg{xk,q)dg{z', Sk) = dg{xk,z'){dg{q, Sk)-dg{xk, q))+dg{xk,q)dg{xk, Sk) 
by the triangle inequality we have 

iT.fTT \i- \ ft \^ frdg{xk,z')dg{sk,Xk)+dg{xk,q)dg{xk,Sk), , , 

^{Uo;xk){xk)~f{sk) < uj{f\— j-, 77— TT ) = ^{f;d{sk,Xk)). 

dg[xk, z') + dg(Xk,q) 

This last inequality clearly implies "^{Uq; Xk){xk) < Uo{xk)- 
Second case, if Sk G V{xk), then 

,T,/rr \/ \ ft \ ^ If dg{xk,Sk)dg{q,Sk) , 

'i'{Uo;xk){xk) ~ f{sk) < sup ^{f;-rY rvTl k)- 

qev(x^) dg{xk,Sk) + dg{xk,q) 



Since 

dg{q,Sk) 



dg{xk,Sk) +dg{xk,q) 



< 1 



On Absolutely Minimizing Lipschitz Extensions 



16 



we have also 'i'{Uo;xk){xk) < Uo{xk)- We conclude that 

Vfc = l,...,7V,^'(C/o;xfe)<C/o . (5.14) 
By Lemma l5. II and this last inequality, we prove inductively that 

Vfc = 1, TV - 1 , ^{Uo; xk+i),..., o ^(L/o; xi) < Uq. 

Therefore we have Ui < Uq and we deduce from Lemma (|5.2ll that sequence 
(J7„)„gN is decreasing. 

By Lemmas 15 . 115 . 2lf01 we prove inductively that 

Un > inf Uois). 

ses 

The sequence {Un)neN is lower bounded and decreasing and therefore converges 
to a function denoted by K{f). It remains to check that H5.1HI . 15.12|I and H5.13|l 
hold. For any e > there exists N ^ N such that 

0< sup iUn{y)~K{f){y))<e ,yn>N . 
yeG-s 

For Xk E G — S and n > N we have 

Un+liXk) = fl{^;Xk) 

with 

* = ^'(C/„; Xk-i) o ^i{U„; Xk-2) o ... o ^'([/„; xi). 
Since ^{y) < C/„(y), Vy £ G - 5 we have 

0<^{y)~K{f){y)<e. 

We can write 

I ^l{K{f)■,Xk)~K{f){xk) \<\ ^l{K{f)■Xk)-^l{^■,Xk) \ + | Ur,+i{xk)-K{f){xk) I . 
Since 

\^l{K{f)■Xk)~ ^x{^■,Xk)\< sup \^{y)-K{f ){y)\<e, 

we obtain 

\^i{Kif);xk)~Kif)ixk) \<2e. 

By lemma we have Uo{s) = /(s), Vs e S* and C/„(s) = /(s), Vs G S, and 
Vn G N. Therefore K{f) is an extension of / and we obtain the stated result. 

□ 
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Now, combining theorems l5.4l and l4.2l functional equation p.l|) has K{f) as unique 
solution. As a consequence of lemmas 15. 115 .215. 31 and of theorems 14. II [4.21 and 

we have the following properties of stability of the extension scheme K: 

Theorem 5.5. Let f , g any two real-valued functions both of domain S . Then 
I K{f){x) - K{f)iy) \< u{f- dg{x, yj) , ^x,y e G ; (5.15) 

sMKU){^) - Kig)ix)) < sup(/(s) - g{s)) ; (5.16) 

xeG s£S 

for any non-empty subsets A, B of S , we have 

sup(X(/ I A){x) - K{f I B){x)) < Mf; Sg{A, B)) , (5.17) 

xeG 

where f \ A and f \ B denote the restrictions of f to A and B and 6g HaussdorfJ 
metric constructed on geodesic metric dg . 

6 Approximation of an AMLE 

Let / denote any f2-continuous real-valued function whose domain is a compact 
non-empty subset S of E. 

In this section we shall consider sequences (G'„,y„)„gN of networks having the 
following properties: 
(Ql) lim„^oo rn = 

where rn ■— sup(i5(G„, E), 6{Sn, S)) and Sn ■— S (1 G„; 
(Q2) lim„^oo Pn = where 

Pn = sup sup d{x,y); 

xeG„ yeV„{x) 

(Q3) lim„^oo II d„ - d 11= 

where d„ denotes geodesic metric on (G„,y„) and 

||d„-d||:= sup \ dn{x,y) - d{x,y) \ . 

x.,y£Gn 

We note bn{x) the open ball of center x E E, radius r„, and Bn{x) the closed ball 
of center x E E, radius p„. 

Lemma l6. II shows that such sequences (G„, Vn)n&i exist in any metrically convex 
metric space. 

Lemma 6.1. ^eguences (G„, V„)„gN exist which satisfy properties (Q1),(Q2),(Q3). 
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Proof. Let (r„)„eN and (pn)neN be any two sequences of positive reals such that: 

Tn <Pn,ne N; (6.1) 

Hm /)„ = 0; (6.2) 
hm — = 0. (6.3) 

Ti->oo pn 

For any x € E, n € N, let us set 5„(x) := {y G E : d{x,y) < r„}, -B„(x) := 
{?/ e : d(a;, y) < pn}- Define (G„, Vn) as follows. Since 5 is a compact subset of 
E, S is covered by balls bn{x), x € S. Therefore there exists xi,...,Xk S S such 
that bn{xi),i = l,...,k cover 5*. Now E — U*L]^6„(a;i) is a compact subset of E. 
Let Xk+i, ■■■,Xm € E - U^^j^bnixi) such that Li'!^^_^_-^^bn{xi) cover E - U*L^i6„(a;i). 
Set Gn ■= {x\, ...,Xm} and, for x G G„, set Vn{x) := G„ fl B„(a;). Note that, 
by construction, wc have S{Sn,S) < Vn and 5{Gn,E) < r„. Therefore properties 
(Ql), (Q2) arc obviously satisfied. 

Now let us show that for n G N sufficiently large we have 

i) {Gn, V^) is a network; 

ii) lim„^oo \\dn- d\\^ 0. 

Properties (PI) and (P2) are immediate. To prove (P4) let x, y € Gn y ^ Vn{x). 
By metrical convexity of E there exists t G E such that d{x, t) = pn — r„ and 
d{x,y) = d{y,t) + d{t,x). Let z € Gn such that d{z,t) < rn- One has d{x,z) < 
d(x,t) + d{t,z) < Pn — Tn + Tn = Pn- Therefore z S Vn{x). Moreover d{y,z) < 
d{y,t) + d{t, z) < Vn + d{x,y) — pn + Tn- Since for n sufficiently large we have 
2T"n — < 0, we infer that d{x, z^ < d{x,y). 

Now let us prove both (P3) and ii). Let iV e N, TV > 1, and x, y G G„. By 
metrical convexity there exists elements of E yo = x,yi, ...,yN = y such that 
d{yi,yi+i) = d{x,y)/N and d{x,y) = Y!^Sq d{yi,yi+i). 

For each i = l, N — 1, choose Zi G G„ such that d{zi, yi) < rn- We have 

d{zi,Zi+i) < d{zi,yi) + d{yi,yi+i) + d{yi+i, Zi+i) < d{x,y)/N + 2rn- 
Now, choosing such that 

d{x,y) 



N 



+ 2r„ < Pn, (6.4) 



we have Zi G Vn{zi+i). Property (PS) is therefore proved. Moreover dg{x,y) < 
Si=o^ d{zi,Zi+i). It follows that 

N-1 N-1 

dg{x,y) -d{x,y) < ^{d{zi,Zi+i) -d{yi,yi+i)) < ^ 2r„ = 27Vr„. 

i=0 j=0 
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Now, for n sufficiently large, one has pn > "irn- Therefore (|6.4|l is satisfied by 
taking = the smaller integer larger than d{x, y)/{pn — "^rn)- It follows that 

2Nr„ < 2d{x, y)/((p„/r„) - 2) + 2r„. 

Therefore |j c?„ — d ||< 2A/((p„/r„) — 2) + 2r„ where A denotes the diameter of 
{E, d). Since lim„^oo Tn — and lim„^oo fyi/ Pn ~ 0, lemma imi is proved. 

□ 

For each n G N, let us define: 

1) the real- valued function /„ of domain S'„ by 

fn{s) = /(s), X e Sn, (6.5) 

2) the real-valued function Wn of domain E by 

W,,{x) = inf {KMn){.s)+uj{f-d{s,x))), x€E; (6.6) 

sSG„ 

where if„(/„) (if„ for short) denotes the solution of (|1.1|1 for network {Gn,Vn) 
under Dirichlct's condition /„. 

Lemma 6.2. We have 

I W„{x) - Wniy) \< u{J-d{x,y)), Wx,y e E ; (6.7) 



"ix^E, inf /(s) < M^„(x) < sup /(s) + r„) ; (6.8) 

s<Edom{f} s£dom{f) 



and 



Knis) - sup c^(/;d„(t,s) -d(t,s)) < Wn{s) < if„(/„)(s), Vs £ G„ ; (6.9) 

teG„ 

/rom which we infer 

Kn{s) - u;{f; \\ d„ - d ||) < M^„(s) < Kn{s), Vs G G„. (6.10) 

Proof. For any x, y ^ E, we have 

W^n(a;) - W^„(2/) < sup (cc>(/; d(s, a;)) - a;(/; d(s, y))). 

sGG„ 

By triangular inequality we have d[s, x) < d{s, y) + d{y, x). 

By growth and subadditivity of uj(f) we have: 

w(/; d{s, x)) < uj{f; d{s, y) + d{y, x)) < w(/; d{s, y)) -f uj{f] d{x, y)). 
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Therefore W^„(x) - W^iy) < Lj{f;d{x,y)). 
For any x ^ E we have 

Wn{x) < infseG„ Kn{fn){s) + sup^gQ^ u}{f; d{x, s)). 

By property of moduU of continuity we have 

Wnjx) < mi s(zG„ Kn{fn){s) + sup^gg /(s) -infsgs/(s)- 

Using (|5.16|) we have inf5gG„ Kn{fn){s) = infsgs„ f{s)- 

Therefore 

Wn{x) < sup^gs /(s) + infsgs„ /(s) - inf^gg /(s) 
and < sup^g^ /(s) + r„). 

Moreover 

W„(x) >inf,gG„i^n(/„)(s). 

Using (|5.16|) again we have infseG„ Kn{fn){s) = infsGS„ f{s)- 
Therefore 

W„(a;) > infsgs /(s). So is proved. 
Let any sq G G„ we have 

t^n(so) - i^„(/n)(so) < /^„(/„)(so) + d(so, Sq)) " i^n(/n)(so) < 0, 

and 

Kn{fn){so) - Wn{so) < SUp^gg^ (ii'™ (/„) (sq ) ^ K„{f„){s) - uj{f;d{sQ,s))). 
By f]-stabilyty of Kn{fn) we have i^„(/„)(s) - Kn{fn){so)) < d„(so, s)) 
Therefore iv:„(/„)(so) - W"„(so) < sup^g^^ (u;(/; d„(so, s)) - uj{f;d{so,s))) 
and 

Kn{fn)iso) - W„(so) < sup cj(/; (i„(so, s) - d(so, s)) < || d„ - d ||). 

So and (IHlTTjl are proved. □ 

Indeed, from Lemma |6.2I sequence {Wn)neN is equicontinuous and equibounded. 
Therefore, by Ascoh's theorem, there exists a subsequence (Wa(n))neN which con- 
verges to a continuous function denoted by u. 

Theorem 6.3. The function u is an AMQ,E of f . 

Proof. We must prove that : 

(i) u is an extension of / 

u{s) = /(s) ,Vs e domif) ■ (6.11) 

(ii) u is r2— optimaly continuous 

I uix) ~ u{v) |< Lo{f- d(x, y)) ,yx, yeE; (6.12) 

(iii) for any open D C E, such that D n dom{f) ^ 0, for any x ^ D, we have 
sup {u{y) ~ uj{u\dD;d{x,y))) < u{x) < inf {u{y) + uj{u \oD;d{x,y))). (6.13) 

yeSD yeSD 
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For typographical convenience let us assume in the proof that subsequence (Wa(n))niEN 
is sequence (W„)nGN itself (the true proof can easily be restated: replace n by a{n) 
almost everywhere). 

Let us show (i). 

For any e > 0, there exists N gN such that Vn > A'', || W„ — u ||oo,b< e- 

For any n > N and s 6 dom(f) there exists by (Ql) Sn G Sn such that d{s, s„) < 

r„. We have 

\uis)-fis) \<Ai+A2+A3, 
where 

Ai -I u{s) - W„(s) I, A2 =\ Wnis) - Kn{s) I, 

^3 =1 K^is) - KniSn) I + I /(S„) - /(S) |. 

We have Ai < e. Using H6.10|l we have ^2 < '^(/; || - ||)- Using we have 

s)) < ^{f\ i"n) + II (i„ — d II). 

In definitive we have 

I u{s) - f{s) \< e + 2Lo{f- r„) + 2^(/; \\ - d ||). 

Since this inequality is true Vn > N and Ve > then, using (Q1),(Q3) and letting 
n tend to 00, we conclude that u{s) — f{s) so we have proved (i). 

The proof of (ii) is immediate by letting n tend to 00 in inequality (|6.7(l . 

Let us show the right inequality of (iii). Let D an open subset of E such that 
D n dom{f) = 0. For any e > 0, there exists N & N such that Vn > iV || W„ - 

U \\oo.E< £• 

Let Ti> N, X e D and y G dD. 
Using (|2.2|l we have 

uj{Wn \aD]d{x, y)) - uj{u \dD\d{x, y)) < 2 || - u ||oo__e- 
Now, setting 

A := u{x) - u{y) - uj{u \gD;d{x, y)), 

we have 

A < 4e + Ai, 

where 

^1 ■.= Wnix)-Wn{y)-L0iWn \aD; d^X , y)) . 

Let £)„ := {z e G„ : 6„(z) n 9^ 0} and (9L»„ = {z e G„ : 6„(z) n ^L* 7^ 0}. 

Since D C UzeD„^n(z) and 9D C Uzga_D„6„(z), we have 6{dD,dDn) < r„ and 

there exists i e -Dn,y G such that 

d(x,x) < S{Gn,E) < r„ and d{y,y) < 6{Gn,E) < r„. 

By lemma we have 

Wn{x) - Kn{i) < u;if; r„) + Lu{f; \\ d„ - d ||) , 

and 

Wn{y) - Kn{y) < coif; r„) + u;{f; || d„ - d ||). 

By (|2.3|l . we have 



c^(M^„ |aD;d(x,2;)) -w(M^„ \9D„;d{x,y)) < MWn;S{dD,dDn)), 
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Since from Lemma [6.21 we have uj{Wn) < ^{f), then 

c^(W„ \9D;d{x,y))^ioiW^ \oD„;dix,y)) <Aioif;S{dD,dD^)) <Aujif;rn). 
Therefore 

Ai<A2 + 6w(/; r„) + 2c^(/; || d„ - d 

where 

A2 := i^„(£) - iir„(y) - uj{Wn \dD„;d{x,y)) . 

Now let us bound 71.2 from above. We write A2 :~ A3 + A4 where 

A3 :=a;(i^„ \9D„;d{x,y)) ~ uj{W„ \gD„;d{x,y)) and 

A4 := Kn{x) - Kn{y) - w(if„ \gD„;d{x,y)). 

We have c?(i, y) — d{x, y) < d{x, x) + d{y, y) < 2r'„. 

Using the subadditivity of a;(W„ |a_D„), we infer that 

LoiWn \gD„\d{x,y)) - Uj{Wn \gD„;d{x,y)) < 2uj{Wn \dD„;rn) < 2w(/;r„). 
Furthermore, using H2.2|l we have 

II \0dJ - i^iWn laoj \\ooM+< 2 || " Wn ||oo,aD„< 2 || if„ - Wn ||oo,G,. • 

Now, using (|6.1U|) we have 

II - W„ ||oo,G„< II - d II). Therefore 

A3<Ai + 2Lu{f; r„) + u;{f; \\ d„ - d ||). 

Now, we bound A4 from above. By theorems 15.41 and 14.21 there exists a unique 
extension v of Kn \dD„ in G„ such that 

( v{z) = fj,{v; z) \fz€Gn-dDn; (rta\ 
\ v{z) = Kn{z) Vz e dD^. ^''■'^> 

Moreover 

- v{q) < uj{Kn \dD„\dn{z,q)), yz,q G G„. 
In particular we have 

v{z) - viq) - Lj{Kn \aD,/,dn{z,q)) < 0, Vg G 9D„, Vz G A.. (6.15) 

Now, we bound sup^^^j^^ | Kn{z) — v{z) \. By symmetry we have only to bound 
A = sxrpz^]j^{Kn{z) — v[z)) frome above. Let 

F = {z G A> : Kn{z) - v{z) = A}, M = supi^„(z), 

zeF 

and 

F = {zGi^:X„(z) = M}. 

Let zq G F be such that 

d{zo, (Sn U dDn)) = inf d(z, (S,, U dD„)). (6.16) 
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Let us first show that we cannot have d{zo, {Sn U dDn)) > 0. Towards a contra- 
diction let us assume it is the case. We have 

v{zq) = ^{v; zq) and Kn{zo) = fJ.{Kn; zq). Using a similar argument to this of the- 
orem ^3 we infer that Vn{zo) C F. Using property (P4) there exists y e Ki(zo) 
such that 

diy, {Sn U dDn)) < d{zo, {Sn U ODn)) 

which is a contradiction with definition (|6.16|) of zq. 

Now, d{zo, {Sn U dDn)) — 0. If Zq G dDn, since Kn{zo) ~ v{zo) we have A = 0. 
If Zq G Sn we remark that zq G Dn- Since zq G Sn H Z)„, Sn r] D = and 

:= {z G G„ : 5„(z) n D 7^ 0} there exists yo & D such that d(zo, 2/o) < 
Moreover since zq ^ D, and by metrical convexity of E there exists po G dD such 
that d{zo,yo) = d(zo,po) + t^(Po,yo) < ^i- 

By definition of dDn there exists go G dDn such that d{po, qo) < ?■„. 
Therefore d{zo, qo) < d{zo,po) + d{po, qo) < 2r„. 
We conclude that 

A = Kn{zo) ~ Kn{qo) + v{qo) - v{zo) < 4uj{f; r„). (6.17) 

From inequalities H6.15|l and H6.17|) . we obtain A4 < 4cj(/;r„). Finally we obtain 

A<Ae+ 12u;{f; r„) + 3Lu{f; \\ dn - d ||). 

Since this inequality is true Vn > N then, using (Q1),(Q2),(Q3) and letting n tend 
to 00, we conclude that 

A < 4e , 

which proves the right inequality of (iii). The proof of the left inequality of (iii) 
is similar but not symmetric because of choice of Wn ■ However it leads to similar 
bounds. □ 



7 Numerical tests. 

The tests of this section are done for the following network: Gn is the set of points 
{i.hjj.h) i,j = 0, ...,n, h — 1/n which densifies SI = [0,1] x [0,1] (eventually 
zoomed and shifted), bn{x) is the ball of center x G Gn radius h, Vn{x) the ball 
of center x G G„ radius k.h. Since norms on are equivalent, balls bn{x) and 
Vn{x) can (and will for convenience of implementation), be choosen to be those 
corresponding to |1 . ||oo or || . ||i. Note that, for fixed n, geodesic metric on 
{Gm Vn) will approach euclidean metric on better and better when k increases. 
The errors in the following tables are 

en,k = sup \ Un,k{x,y) - u{x,y) \ 

where Un,k is the solution of ll.ll of section 1 for S = Sn = dflHGn and f = u \ Sn- 
We first test the algorithm in situations where the solution of the continuous 
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problem is unique and known. u{x, y) = r, 9, r^l'^e^l'^ in polar coordinates, x'*/^ — 
y^/^, in euclidean plane. It is seen that, for a fixed fc, error becomes stationnary 
when n increases. 
Table 7.1: u(x,y) = r 



k/n 


8 


16 


32 


64 


128 


256 


1 


0.023 


0.023 


0.023 


0.023 


0.023 


0.02 


2 


0.0063 


0.0063 


0.0066 


0.0069 


0.007 


0.0067 


3 


0.0062 


0.0031 


0.0031 


0.0031 


0.0032 


0.0032 


4 


0.007 


0.0037 


0.00205 


0.0018 


0.0018 


0.0018 


5 


0.0074 


0.0037 


0.0021 


0.001143 


0.00118 


0.0018 


6 


0.0074 


0.004 


0.0022 


0.001135 


0.000822 


0.00082 


7 


0.0079 


0.004 


0.0023 


0.001178 


0.000602 


0.000571 



(7.1) 



Table 7.2: u(x,y) = 9 



k/n 


8 


16 


32 


64 


128 


256 


1 


0.0251 


0.0141 


0.0139 


0.0138 


0.0138 


0.0138 


2 


0.165 


0.125 


0.0347 


0.00867 


0.00556 


0.00421 


3 


0.236 


0.154 


0.0814 


0.0203 


0.0054 


0.00236 


4 


0.244 


0.191 


0.0958 


0.0347 


0.0088 


0.0012 



(7.2) 



Table 7.3: u{x, y) = r^/^e^/^ 



k/n 


8 


16 


32 


64 


128 


1 


0.156 


0.112 


0.0792 


0.0557 


0.0402 


2 


0.22 


0.159 


0.1123 


0.0812 


0.0557 


3 


0.22 


0.195 


0.1377 


0.0971 


0.068 



(7.3) 



Note that we obtain better approximations if we give "thickness fc/i" to the bound- 
ary that is if we approach the solution of PDE AooU = under Dirichlet's condition 
u \dkHn= uo where de^ = [0, 1] x [0, l]-]e, 1 - e[x]e, 1 - e[. 
Table 7.4: u{x, y) = r 



k/n 


8 


16 


32 


64 


128 


2 


0.0037 


0.0040 


0.0047 


0.0057 


0.0057 


3 


0.0014 


0.0015 


0.0015 


0.0018 


0.0023 


4 


0.0000 


0.0008 


0.0008 


0.0008 


0.00088 


5 


0.0000 


0.0004 


0.0005 


0.0005 


0.0005 



Next we test the algorithm in situations where uniqueness of the solution of the 
continuous problem is not known: Ui{x, y) = x"^ — y^ for || ||oo, 
U2[x,y) =1 x I - I y I for II ||i. 

We note that, in these cases, geodesic metric on Gn coincides, for any k, with metric 
on [0,1] X [0,1]. So, in these cases, we can take k — 1. Numerical tests 17.51 show 
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that the algorithm computes exactely U2 and that error is linear in h for ui (see 
Table 7.5). 

Table 7.5: Ui{x,y) ^ - y^, U2{x,y) =| a; | - | y |. 



n 


8 


16 


32 


64 


128 


256 


el 
e2 


0.06 
0. 


0.03 
0. 


0.015 
0. 


0.0076 
0. 


0.0038 
0. 


0.0019 
0. 



To finish we consider the following two examples. In these examples the metrically 
convex metric space {E,d) is E ^ [0, 1] x [0, l]-]l/4 - e, 3/4 + e[x]l/2 - e, 1/2 + e[ 
for e > small and metric c? on i? is the geodesic metric constructed from local 
euclidean metric. We set uo{x,y) = c?((l/2, 0), (x, y)). Figure 7.1 shows function 
uq in restriction to G„ In the first example we compute the unique solution of 
AooW = 0, u I r = uo where F is the union of the boundaries of internal and 
external rectangles. Numerical tests show that the solution in E is different from 
Uq: this observation corroborates the fact that geodesic cones are not AMLE in 
general metrically convex metric space (see j2|, appendix). 



Figure 7.1: UQ{x,y), n = 100 



'd740.don' 




In the second example we compute A^oU = 0, u \ T — uq where F is now the 
boundary of the external rectangle alone. We obtain a solution which is different 
from Uq and from the solution of the first example. Note the difference between the 
two examples. In the first one we really compute the solution of PDE AqqU = 
because i5 — F is locally euclidean. It is not the case in the second example because 
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the space {E,d) is not locally euclidean at points of the "free internal boundary". 
In fact, in this second example it is likely (we are not insured of the convergence 
of sequence (W^„)„gN in Theorem 6.3) that we compute a AMLE of Uq. 

8 Appendix. 

As announced in Remark 13. 71 we prove that 

11 3 

(- sup u{y) + - inf u(y) - = --Aoou(a;) + o(/i) 

when u is smooth and Du(x) ^ 0. 

Here Ao^u = D'^u[D*u, D*u) and D*u = Du/ \ Du \. 

Since Du{x) ^ then, for h sufficiently small, we have Du{y) ^ for any y G 
Bh{x). Therefore u attains its maximum and its minimum on the boundary 
of Bh{x). Let us denote a;"*" and x~ any points of this boundary such that u{x^) = 
, u{x~) = u~ . Since x"*" is a maximum of u(y) under the constraint \ y — x \— h, 
vectors x'^ — x and Du{x^) have the same direction that is .t+ — x — hD*u{x^). 
For the same reason vectors x^ ^ x and Du{x^) have opposite direction that is 
x~ - X = -hD*u{x~). 

Now, using these expressions of — a; and x^ — x and Taylor formula at x^ and 
x~ we obtain 

2u{x) = u{x+) + u{x^) +A + B + h?o{h) 
where A — h{\ Du \ | Du \ [x^)) and 

B = \h'^{D'^u{x+;D*{x+),D*{x+)) + D'^u{x'' ■D*{x-),D*{x~))). 
Now, using Taylor formula at x for | Du \ {x'^) we have 

I Du I {x+) ^ Du I (x) + hD{\ Du \){x] D*u{x+)) + ho{h). 

By continuity of y ^ D*u{y), it follows that 

I Du I {x+) =1 Du I {x) + hD{\ Du \){x] D*u{x)) + ho{h). 

Similarly, 

I Du I (.T-) =1 Du I (x) - hD{\ Du \){x; D* u{x)) + ho{h). 

Since a straightforward computation shows that D{\ Du \){x] D*u{x)) — Aocu{x), 
and since maps y D*u{y) and y D^u{y) are continuous we obtain in definitive 

2u{x) = it(x+) + u{x^) + 2h'^A^u{x) + /i^Aoou(a;) + h'^o{h), 

which is the announced formula. 
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